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ABSTRACT 

A general two qubit system expressed in terms of the complete set of unit and fifteen 
q traceless, Hermitian Dirac matrices, is shown to exhibit novel features of this system. The well- 
O known physical interpretations associated with the relativistic Dirac equation involving the 

symmetry operations of time-reversal T, charge conjugation C, parity P, and their products are 
?h reinterpreted here by examining their action on the basic Bell states. The transformation 
k-5 properties of the Bell basis states under these symmetry operations also reveal that C is the only 
q\ operator that does not mix the Bell states whereas all others do. In a similar fashion, expressing 
the various logic gates introduced in the subject of quantum computers in terms of the Dirac 
matrices shows for example, that the NOT gate is related to the product of time-reversal and 
^ parity operators. 
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O I. INTRODUCTION 

Qh A system of two qubits A and B offer a canonical paradigm for a large number of 



physical phenomena. The most general physical features of this system have been studied in 



c3 terms of its density matrix tl(A,B), (which is a Hermitian, traceclass, and positive semi-definite 

matrix) written in terms of the parameters of single qubits of A and B and their mutual 
j> correlations. The presence of the mutual correlations in II( A, B) indicate the entanglement 

property of the A and B qubits. Formal results concerning the characterization of the two qubit 
H system have appeared in the literature and a recent account of this may be found in an elegant 
paper of Englert and Metwally [1]. Recently, Rau [2] has constructed propagators for a wide 
variety of qubit interactions. Dorai et al [3] have experimentally demonstrated the usefulness of 
selective NMR pulses in a quantum computation. Arvind and Mukunda [4] have given a two- 
qubit algorithm for finding quantum entanglement. The purpose of the present paper is to employ 
the complete set of sixteen Dirac matrices with a view to delve into possible physical insight that 
may emerge by reinterpreting the existing physical picture that underlies the symmetries of the 
relativistic Dirac theory, namely the time-reversal, charge conjugation, and parity suitably 
reinterpreted for the two qubit system. We follow and adopt the notations given in the review 
article of Good [5] for this purpose. An elegant implication of the algorithm in [4] is found when 
expressed in the Dirac representation and the associated symmetry. 

In section II, we give the mathematical preliminaries and the associated physical 
interpretation of the general 4x4 matrices needed to characterize the two qubit system given in [1 
(see references therein)] by expressing them in terms of the Dirac matrix representation. We give 
in this section an account of the important symmetry operations of time reversal (T), charge 
conjugation (C), and parity (P) and their products introduced in the development of relativistic 
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theory and reinterpret them by working out their actions on the two qubit system expressed in 
terms of the basic Bell states. This leads us to a discussion of the special density matrices given 
in [1] in our language in section III. In section IV, in order to elucidate the symmetry properties 
of T, C, P of the Dirac theory in the context of the two-qubit system, we examine the properties 
of the Bell basis states under the action of these symmetry operations introduced above. These 
reveal several interesting features of the two qubit system. We establish the relationship between 
the Dirac and the Bell representations from which we display the physical significance of the 
Dirac operators in the context of quantum density matrix of two qubit systems. In section V, we 
express some logic gates given for example in [3], in terms of the Dirac matrices and examine 
their properties under the symmetry operations to explore further the underlying physical picture 
that may emerge in this viewpoint. In this section, we also express the results of [4] in the Dirac 
representation. We end the paper with a final section VI, where we make a few concluding 
remarks. 

II. TWO QUBITS IN DIRAC REPRESENTATION - T, C, P 
SYMMETRY 

Each qubit is associated with a set of four 2x2 unit and (traceless) Pauli matrices: 



1 = 
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1 



(o n 
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(\ OA 
~h 
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Two qubits are then associated with sixteen 4x4 matrices formed by direct products of the Pauli 
matrices given above. These are listed in [2]. We note that the Dirac matrices are also formed by 
considering the same direct products. It is also well-known [5] that the Dirac theory of positrons 
and electrons leads to sets of four Hermitian, traceless matrices f^, fl = 1,2,3,4 satisfying the 
relations y y v + J V J^ = 28^ ■ With these four, and their combinations, sixteen 4x4 traceless 

matrices f A such that each of its square is a unit matrix form a complete set in this 4x4 space. 
There are many explicit representations of these but we choose here the representation which 
gives the Schrodinger-Pauli limit of the Dirac theory and connects to Cayley-Klein parameters 
[5]. We list below the equivalent sixteen two qubit matrices with these Dirac matrices: 
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i = i(A)®i(B) 



f, = a 2 (A)®a 1 (B)- 7 2 = a 2 (A)®e 2 (B); y 3 = a 2 (A)® &M> 
f 4 = -a 3 (A)®hB); 



r i(A)®a l (B)^ 



; -'7 3 7i 



J 



r i(A)®a 2 (B)^ 

= S 2 



r i(A)®a 3 (Bf 

= S 3 



if if 4 = °M) ® <^i( 5 ); if if 4 = °M) ® a 2 {B); iy 3 y A = a i (A) ® <7 3 (fi); 

if if 5 = °i( A ) ® if 2 f 5 = ^( A ) ® £ 2 ( 5 ); iy 3 y 5 = ® <^( 5 ); 
/7 4 7 5 = ^(A)®/(5); 



75 = ^(^(8)7(5). 



(2) 



We have arranged the above set to bring out the following known correspondence with the Dirac 
theory: 

the first row is the 4x4 unit matrix; 

the second row of four matrices form a four- vector; 

the third row of six matrices form a second rank antisymmetric tensor; 

the fourth row of four matrices form a four-pseudo-vector; and 

the fifth row is a pseudo-scalar. 

Furthermore, by considering the symmetry properties of the relativistic Dirac equation, the 
following important symmetries were found, which we reinterpret for the two qubit system 
here as follows: 

T - time reversal operator - f 5 f 4 = io 2 {A) ® 1(B); 

C = ch arg e conjugation = -if 2 = C^) ® G 2 {B); and (3a) 

P = parity operator = if 4 = -id 3 (A) ® 1(B); 

Various products of these symmetries are also of interest and they are 

TC = -t 2 ; CP = f 2 f 4 ; Pf = -if 5 ; and 

(3b) 

AAA N ■/ 

TCP = f 2 f 5 . 

In the next section we will use these to construct the density matrices and their transformation 
properties. 
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III. DENSITY MATRIX OF A GENERAL TWO QUBIT SYSTEM 



We begin by representing the 2x2 density matrices of qubits A and B denoted 
respectively by p(A),p(B) with unit trace as follows: 

p(A) = ±(l(A) + ^(A)d ; (A)), p(B) = Ul(B) + Sj(B)&j(B)) 



(4) 



In eq.(4), repeated indices denote the sum over the components, 1,2,3- The 3-vectors with real 
components in (4) parametrize the density matrices. The hat denotes matrix in the above. In 
general these density matrices represent mixed states when the magnitude of the vectors are less 
than unity and pure states when they are unity. In the two qubit system, these density matrices 

are represented by 4x4 matrices with unit trace denoted by 11(A), H(B) and is given by the 
direct products as shown: 



n(A) = i(p(A) ® 1(B)), 11(B) = i(l(A) ® p(B))- 



(5) 



It should be pointed out that these density matrices represent mixed states in the two-qubit space 
even if their parent one-qubit density matrix represents a pure state. In fact, a simple calculation 



with (p 2 (A) = p(A)) in the 2-qubit space leads to (A 2 (A) = -fl(A) . This clearly 



demonstrates the subtle features of the 2-qubit space. 

The 4x4 density matrix formed by a direct product of p(A),p(B) is useful and it is 



Tl d (A,B) = p(A)®p(B) 

r l(A) ® 1(B) + s ; (A)(<7 ; (A) ® 1(B)) + (1(A) ® &j(B)) Sj (B) 
v +s l (A)(a i (A)®& J (B))s J (B) 



(6) 



This direct product density matrix is thus characterized basically by six real numbers plus nine 
others formed out of these as indicated by the last terms in eq.(6), which represent a simple form 
of the correlation between the two qubits in this density matrix. The most general density matrix 
describing two qubits however, has the following structure: 



n(A,B) = 



1(A) ® 1(B) + ^(A)(a ; (A) ® 1(B)) + (1(A) ® ^(B)) Sj (B) 
+ C lJ (A,B)(& l (A)®a j (B)) 



(7) 



Here Qy (A, B) are nine real numbers which are not factorized as in eq.(6), thus representing 

correlation/entanglement between the two qubits. We have chosen the structure of the most 
general two qubit density matrix in the above form so that the marginal density matrices are 
those given in eq.(4): 



4 



Tr B fl(A,B) = p(A) and Tr A fl(A,B) = p(B). 



(8) 



If we subtract eq.(6) from eq.(7) we obtain an interesting expression exhibiting the 
correlation/entanglement between the two qubits explicity: 

fl(A,B) = fl d (A,B) + ^((Qji^B) - Si (Ah(B))(^(A) ® &j(B))) (9) 

We now express eq.(7) in terms of the Dirac matrices by using eq.(2): 



n(A,fi) = 



f l + h (A) y 5 + C 2ll (A, B) y^ + iC^ (A, B) y^ y 5 A 



+s k (B)L k +iC lj (A,B)y J y 4 
where C 24 (A, B) = -s 3 (A), C 34 (A, B) = s 2 (A), and 
K =-iYiYj (i,j,k cyclic). 



(10) 



In this equation the repeated indices are summed with the notation introduced therein. In eq.(10), 
the second term is a pseudo-scalar, third term is a 4-vector, fourth is a pseudo-4- vector, and the 
the last two together form the second rank antisymmetric tensor. Note that our choice of the 
representation of the Dirac matrices explicitly displays the hermiticity of the density matrix. 

It is worthwhile to display the one qubit systems given by eq.(5) in the two qubit space: 

f[(A) = i(p(A) ® 1(B)) = Ul + Sl (A)Y 5 + s 2 {A)iy A y 5 - s 3 (A)h) 



= -(/ + ^(A)/Pf-^ 2 (A)if + s 3 (A)iP), 



(11) 



fl(B) = ±(i(A)®p(B)) = \{i + hB)-Z 

The first expression for the density matrix of the A-qubit involves the pseudo-scalar, and the 
fourth components of the vector and pseudo-vector and in the second equivalent expression for it 
has an interesting structure involving the parity and time-reversal operators. The B-qubit, on the 
other hand, involves the components of the magnetization operator. 
Similarly 



( 



n,(A,5) = 



/ + Si(A)f 5 + is 2 (A)y 4 y 5 - s 3 (A)y 4 + s(B)-Z 
+isM)(hB) • r)r 4 + s 2 (A)(f(B) ■ 7) + is 3 (A)(kB) ■ y)y 5 



(12) 



In the next section, we examine the action of the T, C, P operators on the basic Bell states. 
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IV. T, C, P SYMMETRY PROPERTIES OF BELL BASIS STATES 



The entanglement in the two qubit space is described by the complete set of orthonormal 
Bell-basis states which are given by 




The action of the sixteen Dirac matrices on these states are found to be very useful in elucidating 
the viewpoint we are advocating here and are therefore given in Table I. These four states will 
now be used to rewrite the density matrices considered in the last section from which important 
conclusions concerning concepts of pure and mixed state entanglement will be deduced. By 
using the Bell operator, 

B = 2V2(|3> + )(3> + | - |) = V2(di(A) ® G Y {B) + 03(A) ® & 3 (B)) and its 

square B 2 = 8(|G> + )(<& + | + 1 |) = 4(/(A) ® 1(B) - 2 (A) <g> & 2 {B)\ and 

expressing them in terms of the Dirac representation given in eq.(2), we obtain the following two 
pure state density matrices: 

fl 3 (A,B) = \* + )(*_ 

From the first line in Table I, we deduce the two other pure state density matrices 
ft^A, B) = |ip + )(ip + | = I(/ + iy x y A + f 2 - iy 3 y 5 ), 

1 . (15) 
A 4 (A, B) = I <*>_)(<£_ I = -(/ - iy x y A + y 2 + iy 3 y 5 ). 

The sum of these four is the unit operator in the two qubit space, which is just the completeness 
property of the Bell basis set.Furthermore, twelve other off-diagonal operators of the form, 
I X V+ )( x F m |, I X V+ )(<£+ |, I <&+ )( X F ± |, and \ <&+ )(^> m | are deduced from the transformation 

given in Table I starting from any one of the pure state density matrix given above. This set of 
sixteen density matrix operators based on the Bell states are just another equivalent expression 
for the sixteen operators of the Dirac representation. The relationship of the basic four y- 
matrices with the combinations of the outer products of Bell states is sufficient to establish the 
equivalent relationships between the two and glean from this a physical meaning that can be 
associated with the Dirac matrices in the two qubit context. We therefore list this here: 



I 1 / C . /\ /\ y\ , /\ /\ \ 

I = ^( I - 1 Ti7a- Yi-iWs ). 



(14) 
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- ,T|u/ \/xu 1 




+ 1* VcJ> 


1 - 1* 


V* 11 

/r+ IJ 


Yi 


\\\Ji \/\T/ 1 


1 \T/ \ / \T/ 1 

- T.KT. 




i \rf\ ^ 

+ |4>_ / 


ll 

>(<£_ |J 


r 3 


=4r* + X* + l 


H^-X 4 *- 




+ 4> - 


X*l] 


Ya 






,|+|* + X*. 


_ + \® 


-XM 



(16) 



From these it is interesting to observe that only f 2 has the diagonal Bell density matrices and all 
others represent "decoherent" parts of the density matrix in the Bell basis as they are all of the 
off-diagonal types. Moreover, another important feature is that the four Bell states are the eigen 
states of Y2 ■ 

The above four pure state density matrices represent entangled states because they have the 
following marginal density matrices of A and B qubits, which are maximally mixed chaotic 
states: 

p{A) = U{A)andp{B) = U{B\ (17) 

We also observe that by comparing the four pure Bell-state density matrices in eqs.(14, 15) with 
eq.(10), they have besides the unit matrix, only the correlation terms of the form C u , i — 1,2,3, 
or in the Dirac representation, one part each of 3-vector, pseudo-3-vector, and 2nd rank 
antisymmetric tensor. This is a signature of the presence of an entangled state in a pure state 
density matrix in the two qubit case. 

From Table I, we observe the following important relations: 



1 \ 

Cf¥ ± ) = 


mi|¥ ± ) 


CO> ± ) = ±iO> ± ) 


* 1 , 

P\V ± ) = 




P\* ± ) = -i\* m ) 


f\V ± } = 


±l*m) 


f|* ± ) = ±|*F m ) 


A A. . 


— I * ml 


CP\<S> ± ) = m\<P m ) 


A A, , 

PT\ X ¥ ± ) 


= mi <!>+) 


Pf <£+) = mi ¥+) 


rc|^+X 




7TC|* ± ) = -/|<F m ) 


CPf\ x ¥+ 


H*±> 


y\ y\ y\ . , . 

CPT\<3> ± ) = -\ X ¥ ± 



(18) 



Note from this list, that C is the only operator that does not scramble the states whereas all 
others do. Thus the "charge conjugation" operator seems to have a unique property. 

A third useful representation of the two qubit system is the total spin momentum 
representation, which here consists of singlet with zero spin, and three triplet states with spin 
one with components +1, -1, of its projection along a pre-assigned z-axis. They are explicitly 
represented by 
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s) = \ x ¥_) (singlet); 



In this representation, the pure Bell state density matrices associated with the singlet and the 
triplet with zero projection remain as pure state density matrices and hence without decoherence. 
The other two pure Bell state density matrices have decoherent components in the two projection 
states of the triplet. In the Dirac representation, the unique operator f 2 is decoherent in the two 
projection states of the triplet: 

7 2 = I h X'o I - 1 * X* I - 1 K }(t + 1 - 1 1_ ){t_ |. (20) 

In the following section, we discuss the logic gates etc. in terms of the Dirac basis. 



t n )= 



y* + > + l*-» 



t_) = 



y* + >-|*-» 



(triplet) 



(19) 



V. DIRAC REPRESENTATION OF THE LOGIC GATES IN 
QUANTUM COMPUTATION 



In this section we express several important logic gates in terms of the Dirac 
representation and provide another physical interpretation of these. We also include here a 
discussion of the Arvind-Mukunda algorithm in this frame work which reveals a surprising but 
simple interpretation of their result. 

The NOT and Hadamard yHj gates in the one-qubit representation are 



M)7^ ^ = (^1 + 1^1); 

^■^ + ^)=^[(l T > + WX1 + fl T >-W^I 



(21) 



The two qubit gates are constructed as for example in [2,3], which we re-express in the Dirac 
representation. 



CNOT = XOR = 



(\ 








°1 











1 








1 





v0 


1 








= t(/-i7i7 2 + Ys-iYiYa )• 



(22) 



Note that this contains the operator product of Parity andTime reversal operators. 
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NOT = 



(0 1 0\ 

1 

10 

v 1 



= f 5 = iPT. 



(23) 



This NOT gate is just the product of Parity and Time reversal operators. 



SWAP = 



10 
10 



1 



-(/ + iY 1 y 4 + r 2 + iy 3 Y 5 ). 



(24) 



This operator contains the charge conjugation operator. More interestingly, comparing with 
eqs.(13,14), we find that the SWAP gate is an algebraic sum of the pure Bell states: 



swap = n 1 -n 2 +n 3 + n 4 . (25) 

We now turn our attention to the algorithm suggested in [4] for distinguishing entangled 
from unentangled states, by classifying what they call even and odd density matrices. We first 
express the matrices given by them in our representation. This gives us an insight into the 
algorithm given in [4] . 
We have 



n 



1 



(\ ±i o o\ 

±110 



0) 



-[I-7 A ^iy 2 y,±iy{Y 5 ). 



(26) 



_ 1 



(0 0^ 

1 ±1 

±1 1 





1 / C , /r /r . . yv yv i yv \ 



(27) 



These are both pure state density matrices. We observe that under charge conjugation, tl odd is 

invariant while fl even is not (but fl even is invariant under parity operation but not fl odd ). 
Correspondingly the former is entangled because its marginal density matrix is a mixed state 
while the latter is not because its marginal density matrix remains a pure state. Moreover, by 
comparing with the pure Bell state density matrices (representing entangled pure states) in 
eqs.(14, 15), we find that fl odd is Ajif +sign is chosen, and fl 2 if -sign is chosen. More 
generally, by examining eq.(7), in the light of eq.(27), the signature of entanglement apprears to 
be the presence of only C u and all others absent, as was noted earlier in Sec.IV. On the other 

hand, ft is also a pure state expressed in terms of all the sixteen Bell state density matrices, 
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diagonal and off-diagonal pieces. In fact, we find tl even = | v I / }( v i / |, where 
\ x ¥) = \ x ¥ + ) + \ x ¥_)±\<& + )±\<&_). 

We also observe that the eight even functions encoded by the unitary matrices given in 
[4], in the Dirac representation are 

£[0,4] = / = -£[4,0], U[2,2] = +f 4 , + if 1 f 2 , + if 3 f 5 (28) 

which are all invariant under the parity operation, P. We may also note that they are all separable 
in the sense that they can be respectively expressed as direct products as follows: 

1(A) ® 1(B), (7 3 (A) ®I(B), 1(A) ® d 3 (B), d 3 (A) ® a 3 (B) . (29) 
The odd functions are more complicated and are not invariant under P. This is an alternate way 
to decide whether a given function is even or odd. Thus we have here provided an elegant insight 
into the algorithm given in [4] . 

VI. CONCLUDING REMARKS 

In order to gain clear insight into the physics of any underlying phenomena, the choice of 
a suitable representation is crucial. In the quantum theory of two-qubit systems, the Bell state 
and total spin representations have traditionally been employed with much success in 
understanding the subtle concepts of entanglement etc. In this paper, we have presented one 
more representation in terms of the Dirac Gamma matrices. We then reinterpret the underlying 
symmetries of the Dirac relativistic theory for the quantum two-qubit system. This has elucidated 
and shed light on some recent theoretical works on the subject as is evident from our Sec.V. We 
have also related this to the Bell states and the total angular momentum states to bring out the 
inter-relations among the various representations. Further implications of this representation, it is 
hoped, will be a fruitful avenue of investigation. 
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TABLE I: Action of Dirac Matrices on the Bell states 



The transformation under C, P, T, CP, PT, TC, TCP are shown in parenthesis. 







Yl\ 9 ±)- m \ 9 m) 


y In 1 ) = +\ x ¥ ) 

/2| T ±/ — | T ±/ 


(C = -iy ) 


v ct> \ = Jd) \ 


V hf i = m/lcb \ 

/ 3| ± / m \ ± 1 




y |d> \ = +z'| v T> ) 

l3\^±/ — '| T ±/ 




(£> -x-, \ 

(P = iy 4 j 


/\ 1 \ 1 J. \ 

7 4 |<D ± ) = -|<D m ) 


/\ \\Tt \ 1 1 ^Tn. \ 

7 5 | V P ± ) = ±|$ ± ) 


[PT = -iy 5 ) 


1 JIL \ 1 1 »T( \ 

7 5 |$ ± ) = ±| V P ± ) 


E il*±H*±) 




^ 1 \ 1 \ 

E il*±H*±> 


KK)=-i\^) 


(fc = -t 2 ) 








£ 3l*±H*») 


-ifif4|^ , ±) = H^ , ±> 




- z 'fif 4 *±) = m4> ±) 


-if 2 ? 4 \%) = m\V m ) 


(cp = f 2 f 4 ) 


-^2r4|*±) = ±^ m 






-i7 3 r 4 |*±) = ±W 






W 5 |*±) = | V i , m > 


i?2? 5 \ } ¥±) = -i\®±) 


(cPf = -f 2 f 5 ) 




if 3 f 5 \%)=-\%) 




Z 'f3f 5 |*±) = |*±) 


W 5 | V i / ±) = H*m> 


(f = -f 4 f 5 ) 


W 5 |*±> = H VI/ m) 
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